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We study a model of phantom tethered membranes, em- 
bedded in three-dimensional space, by extensive Monte Carlo 
simulations. The membranes have hexagonal lattice structure 
where each monomer is interacting with six nearest-neighbors 
(NN). Tethering interaction between NN, as well as curvature 
penalty between NN triangles are taken into account. This 
model is new in the sense that NN interactions are taken into 
account by a truncated Lennard- Jones (LJ) potential includ- 
ing both repulsive and attractive parts. 

The main result of our study is that the system undergoes 
a first order crumpling transition from low temperature flat 
phase to high temperature crumpled phase, in contrast with 
early numerical results on models of tethered membranes. 



dered. However, the nature of the crumpling transition is 
still puzzling. Renormalization group (RG) calculations 
H with a Landau continuous model predict a dis- 
continuous phase transition when the dimension d of the 
embedding space is lower than 219, including the physi- 
cal case d = 3. On the other hand, numerical simulations 
of lattice models p|-p^[ , large d expansion and calcu- 
lations based on truncations (20) of the Schwinger-Dyson 
equations are consistent with a continuous phase transi- 
tion. 

In this paper, we try to shed a light on this contra- 
diction with an extensive Monte Carlo (MC) study on 
a model of phantom tethered membranes with bending 
rigidity. As it turns out, our results show that the crum- 
pling transition within our model is of first order in agree- 
ment with the RG prediction. 

Section II is devoted to a description of the model. 
Our method is described in Section III and the results 
are shown in Section IV. Concluding remarks are given 
in the last section. 



I. INTRODUCTION 



Statistical mechanics of membranes is a rich subject 
and has been studied since about twenty years. Motiva- 
tions to obtain a full understanding of the behavior of 
these complex systems are enforced by many experimen- 
tal realizations. For recent reviews, see jl| and [Q]. See 
also H for many introducing and pedagogical courses on 
the subject. 

Membranes are two-dimensional fluctuating systems of 
monomers. According to their physical properties, mem- 
branes can be "fluid" or "tethered" . Fluid membranes 
consist of freely moving monomers, i.e with Hamiltonian 
depending only on the shape of membranes. On the con- 
trary, monomers in tethered membranes are tied together 
by a tethering potential and their connectivity is fixed. 
In addition, a membrane can be self-avoiding if intersec- 
tions with itself are forbidden. Otherwise, it is a phan- 
tom membrane. In this paper, we focus our attention on 
a model of tethered membranes with external curvature 
energy without self- avoidance. 

Any realistic model should include self-avoiding inter- 
actions. But phase diagrams of phantom membranes are 
rich and contribute to understand the behavior of self- 
avoiding membranes It is now firmly established 
that phantom membranes undergo crumpling transition 
between a flat and a crumpled phase. Flat phase pos- 
sesses long-range orientational order between normals to 
the surfaces whereas crumpled phase is totally disor- 



II. THE MODEL 

We consider a 2D lattice of monomers connected in 
hexagonal structure, and embedded in the physical 3d eu- 
clidean space. The tethering potential between nearest- 
neighbor (NN) monomers is a truncated Lennard- Jones 
(LJ) potential. The curvature energy is a standard 
normal- normal interaction between NN triangles. The 
distance between NN monomers are not allowed to be 
larger than an upper bound distance R max - Otherwise, 
in absence of R ma x, monomers are no longer effectively 
tethered at high temperature and the system becomes a 
gaz. In order to keep essential features of the LJ poten- 
tial, Rmax must be sufficiently larger than ro, the NN 
distance corresponding to the minimum of the potential, 
so that Rmax lies in the flat asymptotic region of the po- 
tential. However, to have an actual tethered membrane, 
Rmax should not be too large as discussed above. More- 
over, equilibration times increase as R m ax increases, since 
NN distances are then allowed to grow more and more. 

The system is described by the Hamiltonian 

H= J2 U{r M )-K ]T n a .n p (1) 

<ij> <a/3> 

The first sum is performed on pairs of NN monomers 
< i,j > only, and the second one is restricted to pairs 
of NN triangles < a, (3 >. Tethering interaction between 



1 



IV RESULTS AND ANALYSIS 



2 



NN monomers labeled by i and j depends only on their 
distance r« in the 3d embedding space and is described 
by U(r v ), 



r 7 ) 2 > 



(6) 
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(2) 



with r i; = ||rj — rj||, ri and rj being the position vectors 
in the 3d space. r a is the equilibrium distance between 
NN monomers. The second term in (|l|) is the external 
curvature energy, with K the bending rigidity. The 3d- 
vector iIq, is defined as the normal unit vector of the u—th 
triangle formed by three NN monomers. Note that is 
defined for a counterclockwise oriented triangle. 

The phase space of the model depends on three param- 
eters. We fix two of them, namely (U a and K), and look 
for temperature-dependent properties. 



III. NUMERICAL METHOD 

We consider a membrane of linear size L. The total 
number of monomers is N = L x L. We choose U a = 3 
and K = 1. A more extensive study would require to 
explore the complete phase space, but these particular 
values already give interesting results. The NN distance 
in the ground state is taken to be the unit of distance, 
i.e. r Q = 1, and the upper bound R m ax = 4. We use free 
boundary conditions, at constant pressure in our simula- 
tions. 

The following algorithm was used. Starting from the 
ground state where monomers are on the hexagonal lat- 
tice sites, we heat the system to a temperature T. We 
equilibrate the system at variable volume. The local equi- 
libration is done as follows: we take a monomer and move 
it to a nearby random position in a cubic box of volume 
8 3 around its position, in the 3d space. This position is 
accepted if it lowers the energy. Otherwise it is accepted 
with a probability according to the Metropolis algorithm. 
We repeat this for all monomers: we say we achieve one 
MC step/monomer. We choose 8 ~ 0.1 to have an ac- 
ceptance of the order of 50 %. 

We define the following physical quantities - averaged 
total energy {E), averaged normal vector < n >, aver- 
aged NN distance < d >, radius of gyration R g - with 
the following standard definitions 



(E) = (H) 



(3) 



where < ... > indicates thermal average and the sum in 
< d > is performed only on NN links. 

All the results described below are obtained after ther- 
malization of the system. This requires about 10 6 — 10 7 
MC steps/monomer, depending on the temperature. Af- 
ter thermalization, measures are done on 10 6 -10 7 MC 
steps/monomer, depending also on the temperature. 

Error bars are calculated using a standard jack-knife 
algorithm. 



IV. RESULTS AND ANALYSIS 

At low temperature, equilibrium configurations are 
configurations of minimum energy. This means that NN 
distance must be close to r Q , the minimum of the LJ po- 
tential, and that NN normals must be parallel to mini- 
mize curvature energy. So, equilibrium states correspond 
to a flat state, with < n >~ 1 and R g ~ L. It is an or- 
dered phase. 

At high temperature, maximal entropy configurations 
correspond to crumpled states, where the membrane is 
compact and occupies a very small volume in the em- 
bedding space. It is a completely disordered phase, with 
< n >= and R g < L. 

Between low and high temperatures, a crumpling phase 
transition is expected. In this work, we study the nature 
of this transition to see whether it is continuous or not. 

As a first point, we measured (E), < n >, < d > and 
R g versus temperature for sizes N = 16x16, N = 24x24, 
N = 32 x 32 and N = 48 x 48. Results are shown in 
figures and |. 



C N=16x16 
N=24x24 
□ N=32x32 
* N=48x48 



< n >= 2 (L-1)2 < l ^" a| > (4) 



<d>= j3L-l)(L- 



Wl—T) £ < ry > (5) 



FIG. 1. Averaged energy versus temperature, for 
N = 16 x 16, 24 x 24, 32 x 32 and 48 x 48 
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It turned out that equilibration times are very large for 
this system. So, as we were interested in critical proper- 
ties, we concentrated our work on the temperature region 
around the phase transition, especially for N = 48 x 48. 
Figures p]Mp| and [| clearly show a phase transition be- 
tween a flat and a crumpled phase at T ~ 1.1: there is a 
sharp jump of (E) from the low-T flat phase to the high- 
T crumpled phase. Note that the slope of (E) increases 
with size in the transition region. The order parameter 
< n > drops from a finite value (this value would clearly 
tend to 1 as T tends to zero) to a vanishing value, as 
expected for an order-disorder transition. 



rithmic dependence) in the crumpled phase. It should be 
noticed that the NN distance remains finite in the crum- 
pled phase, which means that monomers actually still 
form a tethered membrane even in the high temperature 
phase as discussed earlier. 
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FIG. 2. < n > versus T, for N = 16 x 16, 24 x 24, 32 x 32 
and 48 x 48 
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FIG. 3. Averaged NN distance < d > versus T, for 
N = 16 x 16, 24x 24, 32 x 32 and 48 x 48 
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FIG. 4. Averaged gyration radius R g versus T, for 
N = 16 x 16, 24 x 24, 32 x 32 and 48 x 48 

The second purpose of our work was to determine 
whether the transition is continuous or not. A standard 
method, when using numerical simulations, consists in a 
finite size scaling analysis of the maximum of the spe- 
cific heat C v . For a second-order phase transition, it is 
expected to grow as the linear system size. For a discon- 
tinuous transition, there are in principle discontinuities 
in thermodynamic quantities. However, for small sys- 
tems, a discontinuous phase transition can appear to be 
continuous if the correlation length is greater than the 
linear size of the system. In that case, C™ ax is expected 
to grow as the size of the system. 

We measured C™ ax for N = 16 X 16, N = 24 x 24 
and N = 32 x 32 using the histogram technique. p2| For 
these small sizes, energy histograms have a single peak 
for all temperatures we explored around the transition 
(multihistogram method). They are found more or less 
gaussian. 

In figure § we plot C™ ax versus N in logarithmic 
scale. Fitting these data, we obtain C™ ax ~ N x with 
x = 0.83(12). This value is far from the value x = 0.5 
expected for a continuous transition. It is closer to 1 , the 
theoretical value for a first-oder transition as discussed 
above. At this stage, in view of this, we conjecture that 
the transition is of first order. As seen below this conjec- 
ture is confirmed by histograms made for a larger size. 



At the transition, R g falls from a finite value dependent 
of the linear size L to a small value more or less indepen- 
dent of L. This corresponds to the scaling R g ~ L v , 
with v = 1 in the flat phase and v = (indeed, a loga- 
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FIG. 5. G™ a * versus N in logarithmic coordinates, for 
N = 16 x 16, N = 24 x 24 and AT = 32 x 32. The dashed line 
is the best fit. 

In order to check further the first-order character of the 
transition, we increase L. For large L, if the transition is 
of first order, the energy histogram should show a struc- 
ture of multiple peaks corresponding to the coexistence of 
ordered and disordered phases at the transition. The sys- 
tem would go back and forth between these phases result- 
ing in a double-peak energy histogram. Taking L = 48, 
we indeed observed this double-peak histogram in the 
region T ~ 1.1, as shown in figure [| This is a very 
strong signal which confirms the first-order character of 
the crumpling transition found earlier by finite-size scal- 
ing of C™ ax . 



V. CONCLUSION 

We have studied the crumpling phase transition by MC 
simulation of a model of tethered membranes with LJ 
potential energy and bending rigidity. We have shown 
clearly the first-order nature of the phase transition be- 
tween flat and crumpled phases, in contrast with earlier 
simulations using gaussian tethering interaction. Note 
that these early MC results |8|-|l^] and also analytical 
calculations ]l9| , p(| show a continuous phase transition 
for models which do not include anharmonic excitations. 
We believe that the anharmonic nature of the LJ poten- 
tial used in our model to some extend contributes to the 
first-order transition observed here. We did not vary in 
this work the value of bending rigidity K. Let us men- 
tion however that for K = 0, there is no flat phase, the 
membrane is crumpled at all T . 

It would be interesting to include self-avoidance be- 
tween non nearest-neighbors in our model. However, it 
seems that in that case the precise form of NN poten- 
tial is irrelevant: for repulsive self- avoidance between non 
nearest-neighbors the membrane is always flat regardless 
of the form of the potential between NN, even in the ab- 
sence of bending rigidity. |p|,p|j2^1 This is interpreted as 
an effective bending rigidity induced by excluded volume 
effect. This can be overcome by an attractive interaction 
between non nearest-neighbors in addition to the repul- 
sive self-avoiding interaction, leading to the folding of 
the membrane at low T and a flat phase at high T. |24| ] 
Including such repulsive and attractive interactions be- 
tween non nearest-neighbors in our model is a formidable 
task which is left for future investigations. 
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FIG. 6. Normalized histogram for the energy, for 
N = 48 x 48. Double-peak structure showing the first-order 
character of the phase transition 
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